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Abstract 

It was recently found that there are classes of nonlocal gravity theories 
that are free of ghosts and singularities in their Newtonian limit [PRL, 
108 (2012), 031101]. In these proceedings, a detailed and pedagogical 
derivation of a main result, the flat space propagator for an arbitrary 
covariant metric theory of gravitation, is presented. The result is applied 
to analyse f(R) models, Gauss-Bonnet theory, Weyl-squared gravity and 
the potentially as}Tnptotically free nonlocal theories. 



1 Introduction 

General Relativity (GR) predicts singularities and doesn't straightforwardly yield to 
quantisation. On the other hand, attempts to modify the theory are restricted by, 
besides phenomenological viability, theoretical consistency. By stra3dng away from 
the Einstein-Hilbert action for the metric of space-time, one easily invites ghosts into 
the theory. While some specific higher derivative theories may be free of ghosts, they 
are not renormalisable, and vice versa [[l[] . 

However, nonlocal theories, featuring an infinite number of derivatives, might pro- 
vide a way around this [Q]. The ultraviolet singularities may then be smoothened 
out without introducing ghosts and while recovering GR predictions at small curva- 
tures. Indeed, the promising attempts at quantum gravity such as string theory and 
loop quantum gravity exhibit nonlocality at some fundamental level. Phenomenolog- 
ically nonlocal gravity has been recently applied to model pre-big bang cosmology 
diiQ], inflation [Jiilal) screening mechanisms [fllsl], dark energy [Igl Eol] . structure 
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formation [IIU Il2f] a nd dark matter nlSi \l4 
renormalisability DlsL and black holes yj 



Theoretical studies have considered 



In these proceedings, we will derive the flat space propagator for the most general 
metric theory of gravity, presented in [Q] ■ In section [2] we write down the action 
and reduce it to a tractable form in the relevant Minkowski limit. Section [3] then 
introduces the formalism and the method to obtain the propagator. The result is 
applied in Section[4]for an analysis of several special cases of interest. 



2 The most General Quadratic Action 



To understand both the asymptotic behavior and the ghost free condition, what is 
relevant is the quadratic action of gravity. In other words if we look at fluctuations 
around the Minkowski background 

then all we need to worry about are the terms that are quadratic in h^^ in the action. 
Now, since in the Minkowski background -R^vAct vanishes, everytime there is a Rie- 
mann tensor in the action it contributes an 0(K) term in the action. Thereby we only 
need to analyse terms in the action that are products of at most two curvature terms, 
and the most general form for the action is given by 

where P is a differential operator containing covariant derivatives and 17^^, . We note 
that if there is a differential operator acting on the left Riemann tensor as well, one 
can always recast that into the above form using integration by parts. 

Since the operator V can only have covariant derivatives and the Minkowski metric, 
one can actually write down the most general action Sq explicitly: 



r 

4 



+ i?^"F5(n)v^v^v,,V;,i?'^^' +J?F6(n)VpV,,V;,v^i?'^^'^'^ +R^;,Fy(n)v,v^R''''^'' 
+ R'''""'''''Pi4iaWpV^v,v^v^v,.VpV^R''''''''] . (3) 
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Using the antisymmetric properties of the Riemann tensor, 

RilJ.v)pcr — -f^/xvCptj) — 0> (4) 

and the Jacobi identity 

one finds after patient index manipulation that the above action reduces to 



So we got rid of 8 of the 14 terms already in a curved background. 



Act 



(6) 



2.1 Linearised Action 

Our next task is to obtain the quadratic (in h^^,) free part of the above action. A very 
important simplification occurs when we realise that the two h-dependent terms must 
come from the two curvature terms present. In other words the covariant derivatives 
must take on the Minkowski values, and we can commute them freely. We then ob- 
serve that the Fg, F13 and F14 terms in the action ([6]) become irrelevant in this limit 
due to the symmetry of the derivative operations contracting the antisymmetric index 
pairs of the Riemann tensor ([4]) . The linearised action contains in the end only 



Furthermore, below it will become clear that there are essentially only two free func- 
tions that determine the properties of the theory in this limit. 

Now our next task is to substitute the linearised expressions of the curvatures in terms 
of h^^: 

1 

V = \i3a^vhl^-d,d^h-nKr^, R=d,d^h>''' -Uh. 

As is obvious, many of the terms simplify and combine. By considering all possible 
consistent contractions, one can deduce that all the terms eventually have to produce 
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an action of the following form: 



1 

2^ 



□ 



(8) 



where we have defined the functions a(n), £>(□), c(n) and d(n) in such a way 
that they reduce in the appropriate limit to the constants a, b, c and d used by van 
Nieuwenhuizen ifigll . The function /(□) appears only in higher order theories. 

We will now compute all of the terms in the original action (O individually. The first 
piece gives 

RF^(n)R = hFiD^h + h^''Fid„d;,d^d„hi''' -hF^nd^d^.h!''' -hi^^F^nd^d^h. 

The third and fourth terms in this case can be combined as follows. Ignoring surface 
terms it is always possible to commute through the local f (□) termfl and we get 

For the two other relevant terms in the ([T]) we obtain in the similar way: 

1 



1 



It remains to relate these terms to the five combinations appearing in action ([8]) . 



2.2 The coefficients in terms of Fi(n) 

In the above section we have calculated the contribution of the higher derivative 
modifications (O to the action ([8]) . We also need to include the contribution from the 
Einstein-Hilbert term, so that the full action we consider is 



S = 



d'^xJ-gR + Sa. 



(9) 



^For non-polynomial terms that is not clear. 
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We then eventually obtain 

1 

a(n) - 1 - -F3(n)n - 2Fio(n)n, 
1 

b(n) - -1 + -F3(n)n + 2f io(n)n , 

1 

c(n)-i + 2Fi(n)n + -F3(n)n, 

d(n) - -1 - 2fi(n)n - ^FgCn)^, 

/(□) - -^Fi(n)n - ^F3(n)n - ^Fio(n)n . (lO) 

From the above expressions we observe the following interesting relation^ 

a + b = 0, c + d = 0, b + c+f=0, (11) 

so that we are really left with two independent arbitrary functions. This can be un- 
derstood as a consequence of the Bianchi identities, as will be shortly clarified. 



2.3 Field Equations & Bianchi identities 

The field equations can be derived straightforwardly by varying the action ([8]) : 

+r,^,d(n)nh + f[n)n-'d^d^d^dy'' = -KT^,. (12) 

The matter side is conserved by the stress energy conservation and the geometric part 
because of the generalised Bianchi identities [20] due to diffeomorphism invariance. 
Thus 

- /ctV^t^ = - (c + d)nd,h + (a + b)nh':.^ + ib + c+nh'^f^^ . (13) 

It is then clear why (fTTl) had to hold. 

The above field equations can be written in the form 



where F^^J'^'^ is the inverse propagator To compute the propagator, we need to learn 
to deal with spin projector operators. 



^An immediate observation we can make is that one cannot construct the Fierz-Pauli term, for which 
a = — d ~ and b = c = from an action hke Massive gravity is not among the metric theories 
we consider here, but inherently bimetric. 
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3 Propagators 

In this section we are going to derive the propagators for the field equations (I12D . The 
basic algorithm is as follows [fl9[l : First, we express the field equations in the form 
gi4D . where the inverse propagator is expressed in terms of the six operators, Vi'. 

6 

n-^ = Y,CiVi (15) 

i=l 

In the momentum space description of the coefficients Cj's are scalars which can only 
depend on k^. Finding the suitable operators, it is possible to decompose the field 
equations into a decoupled set of equations of motion for the relevant degrees of 
freedom. These are then conveniently invertible. 



3.1 Spin projector operators 

Let us introduce 

2 _ 1 1 

^ _ 1 

^ — l^^np^va-^^is.a^vp-^^vp^tia^^va^tip^^ 

= ^^^Mv^pa. Vl=^o,^yQ,,, (16) 

where the transversal and longitudinal projectors in the momentum space are respec- 
tively 



Note that the operators are in fact 4-rank tensors, 'P^^p^., but we have suppressed 
the index notation here. 

Out of the six operators four of them, {V'^,V^,V^,V^}, form a complete set of pro- 
jection operators: 

K'Pi = 5''5abK and V^ + V^+V° + V°^l, (17) 

as one can easily verify. These projection operators together represent the six field 
degrees of freedom, the additional four fields in a symmetric tensor field, as usual, 
being gauge degrees of freedom. V'^ and represent transverse and traceless spin 
2 and spin-1 degrees, accounting for four field degrees of freedom, while V^, 
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represent the spin-0 scalar multiplets. In addition to the above four spin operators we 
have V^^ and V^^ which can potentially mix the two scalar multiplets. In particular, 
we have that 

-Pfj-Pk-Sj.V^j' T^'K-SaS^.Vl n°^°=5,,p;., (18) 
as one may again easily check. 

From (fT5]) and (ITtI ) it trivially follows that we can write ([141 ) as 

6 

Y,CiVih = KliV^ + V^+V° + Vl)T. (19) 

i=l 

By multiplying with the different projector operators on either side of the equation 
we can now obtain the decoupled field equations for the different spin multiplets. 

3.2 Inverting the field equations 

Having outlined the algorithm for finding the propagators, let us now proceed to 
obtain them in our model specified by the action ([8]). We need to express all the 
operators in (I12D in terms of the operators . Let us start with 

One can continue in an analogous fashion to obtain the projector decomposition of all 
the operators appearing in the field equations (I12D . For the first three terms we then 
obtain 

b(n)a,%,h- ^ _b(_fc2)fc2[pi + 2P°]h, 

While all the above operators appear in two derivative generalisations of gravity and 
were discussed in [|19f] . the / term in (I12D is specific to higher derivative theories. Its 
decomposition is rather simple 

We are now ready to write down the projected field equations, and the corresponding 
propagators. By acting with "P^ on ( |19D we find 

ak^V^h = KV^T^V^h = K\ ^^T. (20) 
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Similarly, acting with V^, one finds 

(a + b)fc2pi/, = Kp^T. (21) 

Rather interestingly, since recalling Eq. (IllD we know that a + b = 0, this implies 
that there are in fact no vector degrees of freedom, and accordingly the stress-energy 
tensor must have no vectorial part: V^t = 0. 

Next let us look at the scalar multiplets. By acting and on gi9D we obtain 

ia + 3d)k^V^}i + {c + d)k^V3V°^h = kV^t, (22) 
{c + d)k^V3r°^h + ia + 2b + 2c + d+f)k^Vy = kV^t . (23) 

As we see, in principle, the scalar multiplets are coupled. However, by applying the 
projector on equation (I22D or the projector on equation (I23D from the right 
hand side, one sees that c + d = 0, in accordance with (IllD . The scalars decouple and 
one can now straightforwardly invert the field equations to obtain the propagators: 



(a + 3d)k^r°h = kV\ ^ V% - k- ^— and (24) 



s 



(a + 2i, + 2c + d + 4/))c^P> = ^ Vy = .^-^^-^-^-j^^, . (25) 

respectively. The denominator corresponding to the P° projector vanishes. So there 
is no w-multiplet, but the 5-multiplet picks up a nontrivial propagator. 

To finally summarise: 

"=5?+rap- '''' 

We have thus arrived at the main result of [Q] ■ 
4 Applications to special cases 

In this section, we consider the implications of the result (1261) to some special cases. 



4. 1 General Relativity 

Since we want to recover GR in the infrared, we require from any viable theory that 

a(0) = c(0) = -b(0)--d(0)-l, (27) 
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corresponding to the GR values. In GR these functions are the same constants for any 
Fourier mode. The above condition ensures that as — > 0, we have only the physical 
graviton propagator, 

lim n = iV^k^) - (P°/2/c2) = Hgr . (28) 

There is a crucial subtlety one should observe here. Although the residue at fc^ = 
is negative, that is a benign ghost. In fact, has precisely the right coefficient to 



cancel the unphysical longitudinal degrees of freedom in the spin-2 part []19f1 



4.2 Gauss-Bonnet gravity 

Let us consider the theory C = R + a(n)G, w^here G is the Gauss-Bonnet invariant 
G = - 4R^i,R'^'' +R^yp^R'^''P^ and the function a(n) in the simplest case can be 
just a constant coefficient. As is well known, in four dimensions the Gauss-Bonnet 
term is a topological invariant that does not contribute to the gravitational field equa- 
tions. Therefore it is not a surprise that it doesn't introduce any modifications to the 
propagator either In Eq.© we have now f ^ = a, F3 = —4a and F^q — ^t- Regardless 
of a we then see from Eqs. JlOD that a = c = — fa = —d = 1, and thus the properties of 
the theory are identical to GR. 



4.3 L(R) gravity 

The C(R) gravity is a popular subject of study. For our purposes here, it is enough to 
consider the expansion of the lagrangian around flat space, 

OiR) = £(0) + C'iO)R + ^C"iO)R^ + ■■■ . (29) 

The first term one identifies with the cosmological constant, £(0) = —2k~^A, and the 
second term should reduce to the Einstein-Hilbert term in a viable theory, C'iO) = 1. 
The relevant modification of the theory is then given by the quadratic part. Since 
only Fi is now nonzero in ([3]), we readily see from (IIOD that then a = —b = 1, 
c = -d = 1 - C"{_0)n and / = -C"{,0)n. The propagator is thus 

-p2 po 

n — ^ (30) 

k^ k\2 + 3£"i0)k^)' 

The scalar part of the propagator is modified. Since these theories are a specific class 
of scalar-tensor theories, we expect an extra scalar degree of freedom. Its appearance 
can be made transparent by rewriting the above result as 
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Indeed, the C{R) correction entails an additional spin-0 particle which is nontachyonic 
as long a^ £^^(0) > 0. One also sees that though these theories are classically viable, 
they cannot improve the ultraviolet properties of GR since the graviton propagator 
retains its form. 



4.4 Conformally invariant gravity 

As an example of a ghastly theory, let us consider the Weyl squared gravity. The Weyl 
tensor is defined as 

R 1 

The theory is then specified by the conformally invariant Weyl-squared term, C = 
R — Kc^, w^here m is the mass scale at which the correction becomes relevant. It is 
Straightforward to compute that 



- Vpa^''"''" - ^^vR"" + , (32) 



from which we quickly infer, using again Eqs. (l3llOD that now a = —b = 1 — (k/m)^, 
c = -d = 1 - [klrnf 13 and / = -2(fc/m)^/3. We obtain the propagator with a 
double pole for the graviton: 



-p2 -pO -p2 

F (1 - (fc/m)2) ~ 2fc2 = ~ • (33) 



From the latter form of the propagator it is obvious that the theory contains an extra 
spin-2 degree of freedom with respect to GR. Moreover, the extra contribution always 
comes with the wrong sign: this is the Weyl ghos10. 



4.5 As5Tnptotically free gravity 

Finally, we show how the ultraviolet properties of GR are improved via a nonlocal 
extension. Just for simplicity, let us restrict to the special class of theories with f = 0. 
From (IllD we then see that a = c = —b = —d. This means that we are essentially left 
with just a single free function 

a(n)-i-iF3(n)n-2fio(n)n- (34) 

^For an alternative derivation and generalisation of this stability condition, see [|2l|] . 

"•However perhaps the negative norm states can be consistently projected out of the Hilbert space [|23] . 



10 



Barcelona Postgrad Encounters on Fundamental Physics — Trim Size: 10" x 7" — "entry-koivisto.tex" 



Prepared for "Proceedings of the Barcelona Postgrad Encounters on Fundamental Physics" — 2013/2/5 — #11/?? 



Thus, the GR propagator is now modulated by the /c-dependent function a(n). We 
now reahse that as long as a(n) has no zeroes, these theories contain no new states 
as compared to GR, and only modify the physical graviton propagator. Polynomial 
functions would correspond to higher - but finite - order gravity, and would inevitably 
result in new pathological states. This can be avoided in nonlocal, i.e. infinite or- 
der higher derivative theories. Furthermore, by choosing a(?c^) to be a suitable entire 
function we can indeed tame the behavior of the ultraviolet gravitons. A simple exam- 
ple can be provided by a = exp (fc/M)^, where M is a mass scale at which the nonlocal 
modifications become important. The integrals over the propagator quickly tend to 
zero at high momenta k> M and we expect finite results from physical calculations. 

5 Conclusions 

Having derived the main result (1261) . we considered its implications in some special 
cases. We readily reproduced the known results: while GR and Gauss-Bonnet theory 
share the same field content, / (i?) gravity has an extra healthy scalar and Weyl gravity 
an extra pathological spin-2 field. New classes of nonlocal theories were found, that 
are both unitary and devoid of singularities. The ongoing further work includes the 
generalisation of the result to curved backgrounds. 
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We obtain a very simple expression for the propagator: 




(35) 
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